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Abstract 


This  thesis  concerns  the  development  of  a  finite 
element  model  of  the  T-38  horizontal  stabilator  for  use 
on  NASTRAN.  The  model  is  to  be  used  to  analyse  degrada¬ 
tions  in  flutter  speed  due  to  repair. 

Static  analysis  has  shown  the  model  to  be  lacking  in 
torsional  stiffness.  The  probable  cause  being  the  inability 
of  NASTRAR  plate  bending  elements  to  model  torsion  .'.  '11s. 

An  Increase  of  elastic  and  shear  moduli  of  plate  bending 
elements  in  the  model  by  30  percent  produced  more  accurate 
results  but  additional  investigation  is  necessary. 

Modal  analysis  has  pointed  to  a  modeling  error  in  the 
root,  trailing  edge  area.  The  affect  has  caused  an  addi¬ 
tional  node  to  appear  on  the  trailing  edge  for  modes 
above  1.00  cps  in  a  free-free  condition.  Investigation  of 
the  steady  aerodynamic  pressure  distribution  over  the 
stabilator  shows  good  correlation  with  experimental  results. 

A  flutter  analysis  procedure  was  established  and  the 
affects  of  the  errors  found  in  the  structural  model  were 
investigated.  With  no  corrections  made  to  the  model,  a 
flutter  speed  equivalent  to  that  predicted  using  strip  theory 
was  achieved  for  the  sea  level  condition. 
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INVESTIGATION  OF  AN  IMPROVED 


FLUTTER  SPEED  PREDICTION 
TECHNIQUE  FOR  DAMAGED  T-38 
HORIZONTAL  STABILATORS  USING  NASTRAN 


INTRODUCTION 


Throughout  the  life  of  the  T-38  the  area  subject  to  the 
most  scrutiny,  from  an  aeroelastic  viewpoint,  has  been  the 
horizontal  stabilizer /elevator  (stabilator  or  just  stab). 
Several  problems  exist  which  effect  the  mass  distribution  of 
the  airfoil  and  change  its  elastic  characteristics.  Most 
prominent  are  the  effects  of  water  absorption,  skin  delamin¬ 
ation,  and  cracks.  Resulting  changes  in  mass  distribution, 
directly  or  as  a  result  of  repair,  effectively  degrade  the 
flutter  velocity. 

To  date,  repair  limits  have  been  based  on  calculations 
made  using  Strip  Theory  (Ref  1,2).  In  the  two  years  prior  to 
July  1980,  sixty-six  port  side  and  thirty-six  starboard  stab- 
ilators  were  purchased  as  replacements  at  the  cost  of  approx¬ 
imately  $10,000  each.  Only  three  stabilators  in  that  period 
were  repairable  (Ref  3).  It  is  hoped  that  a  finite  element 
model  can  be  used  for  more  accurate,  less  conservative  predic¬ 
tions  and  that  more  of  the  stabs  can  be  saved.  Probable 
savings  will  be  proportional  to  the  increase  in  accuracy 
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achieved  with  the  new  model. 

This  project  is  a  follow  on  to  the  thesis  by  John  0. 
Lassiter,  AFIT/GAE/AA/80M-2  (Ref  4),  who  successfully  con¬ 
structed  and  partially  verified  a  finite  element  model  of  the 
stab  for  use  on  NASTRAN.  Several  problem  areas  were  identi¬ 
fied  by  Lassiter  and  are  topics  of  investigation  in  this  thesis. 

Two  designs  of  the  stab  were  proposed  in  the  initial  de¬ 
velopment  of  the  T-38,  referred  to  here  as  the  Series  2  and 
Series  3  models.  Experimental  information  was  compiled  on  the 
Series  2  model  only,  while  the  Series  3  model  was  the  chosen 
design.  For  this  reason,  deflections  of  a  model  of  the  Series  2 
stab  under  static  loads  is  the  topic  of  Chapter  1.  Other 
chapters  include  investigation  of  aspects  critical  to  accurate 
flutter  speed  prediction  for  the  Series  3  stab.  Chapter  2 
involves  a  vibration  analysis  of  the  Series  3  stab  both  exper¬ 
imentally  and  analytically.  Some  discussion  of  the  vibration 
characteristics  of  the  Series  2  stab  is  also  included  to  fully 
integrate  the  relationship  between  the  two  models.  Chapters 
3  and  4  depart  from  structural  aspects  to  consider  steady  flow 
over  the  airfoil  and  the  importance  of  wing/ fuselage  effects  in 
the  flutter  analysis  respectively. 
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CHAPTER  1 


VERIFICATION  OF  THE  SERIES  3  FINITE  ELEMENT 
MODEL  UNDER  STATIC  LOADS 


1.1  Introduction 


Discrepancies  between  theoretical  results  from  the 
Series  3  finite  element  model  and  test  results  reported  in 
Ref  4  suggest  further  investigation  and  verification  of  the 
structural  and  aerodynamic  models  of  the  Series  3  T-38  stabi- 
lator  (Fig  1.1).  Because  test  information  has  been  compiled 
for  the  Series  2  stab  only,  these  errors  can  be  attributed  to 
two  sources.  Either  the  Series  3  stab  has  been  modeled  in¬ 
correctly  and/or  there  is  an  appreciable  difference  between 
the  Series  2  (Fig  1.2)  and  Series  3  design. 

It  was  assumed  in  Ref  4  that  the  finite  element  model  of 
the  Series  3  stab  was  in  error  because  of  a  statement  made  in 
NAI-57-59  (Ref  5)  that  the  Series  2  and  Series  3  stabs  have 
essentially  the  same  stiffness.  To  investigate  the  accuracy 
of  this  statement,  a  model  of  the  Series  2  stab  was  constructed. 

Because  resonant  frequency  and  flutter  speed  are  both 
directly  related  to  stiffness,  a  great  deal  of  attention  was 
focused  on  results  that  suggested  the  stiffness  was  incorrect. 
Experimental  results  from  an  influence  coefficient  study  and 
static  load  test  of  the  Series  2  stab  (Ref  6)  were  the  main 
source  of  comparison.  All  numerical  analysis  for  this  section 
was  accomplished  using  NASTRAN  Rigid  Format  1  (Ref  7). 
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1.2  Static  Load  Comparison  of  the  Series  2  and  Series  3 


Stabilators 


1.2.1  Analytical  Models 


Figures  1.1  and  1.2  show  the  most  prominent  differences 
between  the  Series  2  and  Series  3  stabs.  Reference  4, 

Appendix  A  discusses  the  FORTRAN  program  BDATA  that  was  cre¬ 
ated  and  used  to  generate  the  NASTRAN  bulk  data  deck  for  the 
Series  3  model.  For  an  extensive  discussion  of  the  Series  3 
construction  and  finite  element  model  see  Ref  4:8-14.  Addi¬ 
tional  modules  added  to  BDATA  include  generation  of  'CBAR' 
and  corresponding  ' PBAR '  cards  that  represent: 

1.  the  leading  edge  extrusion 

2.  the  tip  rib 

3.  the  torque  tube 

4.  the  auxiliary  spar  * 

5.  and  the  intermediate  ribs  * 

Physical  dimensions  used  to  construct  these  modules  were  taken 
from  blueprints  reproduced  from  aperture  cards  (Ref  8). 

Two  models  of  the  Series  2  stab  were  generated  using 
BDATA.  One  of  the  models  was  used  in  an  eigenvalue  analysis 
and  will  be  discussed  later.  The  other  is  an  influence  co¬ 
efficient  model  that  contains  additional  grid  points  that 

*  These  items  are  included  only  In  the  Series  2  stab 
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Figure  1.2,  Series  2  Stabilator 


correspond  to  points  of  load  applications  and  observed 
responses.  Figure  1.3  is  an  exploded  view  of  the  Series  3 
model.  Quadrilaterals  represent  NASTRAN  ' CQUAD1 '  elements 
while  dashed  lines  represent  ' CBAR*  elements.  figure  1.4 
shows  the  node  numbering  convention  for  all  the  models. 
Figures  1.5  and  1.6  show  the  two  Series  2  models  used.  Tri¬ 
angles  in  the  Series  2  influence  coefficient  model  represent 
NASTRAN  ' CTRIA1 '  elements. 

1.2.2  Loads  and  Boundary  Conditions  for  the  Series  2 
Influence  Coefficient  Model 

Reference  4:24-28  discusses  the  single  point  constraints 
which  were  applied  to  the  model  to  simulate  conditions  of 
symmetry.  These  constraints  were: 


GRID  POINT 

DISPLACEMENT 

ROTATION 

140 

145 

X  Z 

Y-AXIS 

146 

Y 

X , Z-AXIS 

Constraining  grid  point  140  about  the  Y-axis  removes  any  rota¬ 
tion  that  may  have  been  induced  through  the  actuator  assembly 
or  additional  rotation  from  torque  tube  twist  (Ref  4).  Con¬ 
straints  on  grid  point  145  are  typical  of  a  bearing  and  those  on 
146  are  used  to  invoke  symmetry.  Table  I  presents  the  static 
loads  and  points  of  application. 
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Lgure  1.3,  Series  3  Model,  Exploded  View 
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TABLE  I 


Influence  Coefficient  Study,  Loads  and  Loading  Points 


LOAD  (lb) 

HSS 

7o  CHORD 

FIGURE 

400 

82.0 

52.7 

1.7 

400 

70.75 

20.0 

1.8 

400 

70.75 

75.0 

1.9 

1200 

42.25 

20.0 

1.10 

1200 

42.25 

75.0 

1.11 
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1.3  Results  and  Conclusions 


Figures  1.7-1.11  show  a  comparison  of  the  deflections  of 
the  Series  2  stab  and  Series  3  stab  under  identical  static 
loads  and  boundary  conditions.  In  these  figures,  solid  lines 
are  NASTRAN  results.  Dashed  lines  are  NOR-60-6  results. 

Blacked  in  symbols  and  open  symbols  represent  Series  2  and 
Series  3  results  respectively. 

It  is  easily  observed  that  the  Series  2  model  approximates 
the  experimental  results  more  closely  than  the  Series  3  model 
but  is  still  lacking  somewhat  in  stiffness.  Remember  that  the 
experimental  results  were  obtained  for  a  Series  2  and  not  a 
Series  3  stab. 

These  results  suggest  a  lack  of  stiffness  throughout  the 
model.  In  an  attempt  to  make  an  overall  correction,  the  elas¬ 
tic  and  shear  moduli  of  the  material  comprising  the  majority 
of  the  stab  were  increased  by  30  percent.  The  figure  30  per¬ 
cent  will  be  discussed  later.  Parts  of  the  structure  affected 
were : 

1.  the  skin 

2.  outboard  30  percent  of  the  main  spar 

3.  the  auxiliary  spar 

4.  and  the  intermediate  ribs 

Noticeably  absent  from  the  list  are  the  torque  tube,  the 
majority  of  the  main  spar,  and  the  root  rib.  These  provide 
spanwise  bending  stiffness.  The  elements  listed  affect  an 


13 


o 

o 


in 

«*■» 


\  ™ 
\  °° 
\  w 
\  = 


,  to 

\  J2 


9 

3 

O 

r*0 

CN 

vD 

u 

O 

_c 

tn 

in 

O 

4-> 

O 

<y  vX) 

•— 

•— 

1 

in 

Lmm 

■ 

u 

v. 

c£ 

<D 

*1 

1 

<u 

O 

O 

•— 

CO 

CO 

Z 

• — 

a. 

1 

1 

1 

CL 

d) 

m 

<u 

CD 

in 

*•— 

in 

•o 

CD 

"O 

“D 

O 

0) 

c 

“O 

L. 

LU 

0 

.a 

c 

0 

Ui 

T3 

0  *0 

X) 

5 

•— > 

CD 

u 

.n  u 

a> 

e 

> — 1 

O 

cn 

0 

0  O 

c 

>*  CO 

• — 

c 

-C 

-C 

•— 

GO 

“O 

.— 

etf  0 

•— 

Z) 

in 

“O 

— 

c 

0 

X 

•— 

<0 

•  cN> 

0 

a 

<u 

IS) 

-C 

<u 

O 

C\i  lf\ 

L. 

Q. 

«— 

0 

h- 

CNl 

i/NSh 

O 

CD 

0 

x— ' 

(saqoui)  uoT^ooxjaa 


at  HSS  30.00 


Percent  Chord  at  HSS  30.00 


increase  in  chordwise  bending  stiffness.  Appendix  A  presents 
a  comparison  of  the  Series  2  and  Series  2  stiffened  models 
with  the  experimental  results.  Figures  A1-A5  show  the  stif¬ 
fened  model  has  approximated  the  experimental  results  much 
more  closely  in  pure  bending.  The  largest  errors  result  from 
loads  applied  off  the  elastic  axis  (torsion  loads).  Further 
'tuning'  of  the  model  is  necessary  to  alleviate  this. 

It  can  be  concluded  that  comparing  the  displacements  of 
the  Series  3  model  with  experimental  data  for  the  Series  2 
stab,  as  was  done  in  Ref  4,  is  not  an  accurate  approach.  It 
can  also  be  concluded  that  both  the  Series  2  and  Series  3 
models  are  lacking  in  overall  stiffness  probably  more  chord- 
wise  than  spanwise.  It  is  probable  that  the  error  arises  from 
the  inability  of  NASTRAN  flat  plate  elements  to  model  torsion 
cells.  Remember  the  model  is  composed  solely  of  bar  elements 
and  flat  plates.  In  view  of  these  findings  the  plate  elements 
may  be  given  additional  stiffness  to  compensate.  The  next 
chapter  contains  evidence  which  justifies  this  further. 


CHAPTER  2 


MODAL  ANALYSIS  OF  THE  SERIES  3  STABILATOR 


2 . 1  Introduction 

Throughout  this  chapter ,  any  mention  of  a  stabilator  will 
refer  directly  to  the  Series  3  stab  unless  otherwise  stated. 

To  further  verify  the  finite  element  model,  a  comparison 
of  natural  frequencies  and  mode  shapes  was  desired.  Three 
experimental  methods  for  determining  frequency  response  infor¬ 
mation  are  discussed  in  the  following  sections.  Each  method 
is  introduced  followed  by  some  discussion  of  the  test  setup, 
data  acquisition,  and  results. 

In  each  test  the  stab  was  supported  by  elastic  chords, 
attached  at  each  corner  with  adhesive,  and  hung  parallel  to 
the  floor  (see  Fig  2.1).  Constraints  oE  this  type  represent 
a  Free-Free  boundary  condition.  Table  II  is  a  list  of  the 
equipment  used  in  each  of  the  tests.  Any  mention  of  a  force 
gauge,  accelerometer,  etc.  in  the  following  sections  refers 
directly  to  those  items  listed. 
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Figure  2.1,  Stabilator  Vibration  Test  Setup 


2.2  Sine  Dwell  Test 


2.2.1  Test  Setup  and  Equipment 


A  shaker  was  attached  to  the  underside  of  the  stab  at 
Fuselage  Station  (FS)  496.25  and  Horizontal  Stab  Station 
(HSS)35.25.  A  force  gauge  was  located  between  the  shaker  and 
stab  and  an  accelerometer  was  placed  at  the  outboard  tip, 
trailing  edge.  Outputs  from  the  force  gauge  and  accelerometer 
were  amplified  and  passed  to  a  dual  trace  oscilloscope  for  ob¬ 
servation  of  phase  and  variation  in  amplitude.  Output  from 
the  accelerometer  also  ran  through  a  frequency  counter. 


TABLE  II 


Experimental  Apparatus 


APPARATUS 

MANUFACTURER 

MODEL 

COMMENTS 

Fourier  Spec¬ 
trum  Analyzer 

Hewlett- 

Packard 

HP  5451B 

25  lb.  shaker 

Ling 

Electronics 

Accelerom¬ 

eters 

Vibrametrics 

1000  A 

Force  Gauge 

Vibrametrics 

208  A03 

1  mV/lb  force 

Amplifiers 

Intech 

A2318 

Variable  gain 

Oscilloscope 

Hewlett- 

Packard 

HP  1707B 

Universal 

Filter 

General- 

Radio 

Bandpass 
50-1000  Hz 

Terminal 

Tektronix 

TEK  4010-1 

Copier 

Xerox 

Versatek 

Hard  copies 
from  terminal 

Voltmeter 

NLK 

LX- 2 

, 

Power  Supply 

AFFDL/FBG 

±  15V  DC 

Force  Gauge 
Power  Unit 

Piezotronics 

480  A 

Frequency 

Counter 

Hewlett- 

Packard 

HP  5245L 

2.2.2  Data  Acquis ition  and  Results 


Frequencies  were  observed  to  be  resonant  when  large  in¬ 
creases  in  the  amplitude  from  the  response  accelerometer 
occurred.  Each  resonant  frequency  encountered  was  also  easily 
audible.  When  a  resonant  frequency  was  observed  it  was  recor¬ 
ded  and  an  attempt  was  made  to  identify  the  mode  of  vibration 
(for  lower  frequencies)  by  touching  the  surface  and  looking 
for  node  lines.  Table  III  contains  the  resonant  frequencies 
found  from  0-300  cps. 

Resonant  frequencies  with  node  lines  over  the  shaker 
attachment  point  could  not  be  identified.  These  were  deter¬ 
mined  using  the  other  methods  presented  in  sections  2.3  and 

2.4. 

An  advantage  to  the  Sine  Dwell  method  was  that  all  reson¬ 
ant  modes  found  were  audible  while  some  were  not  easily  detec¬ 
table  on  the  oscilliscope .  The  audibility  factor  became  use¬ 
ful  when  line  noise  at  60  cps  and  180  cps  (due  to  the  use  of 
3-phase  power  in  the  lab)  made  electronic  detection  of  reson¬ 
ance  near  these  frequencies  impossible. 

2 . 3  Modal  Assurance  Criteria 

2.3.1  The  Modal  Assurance  Criteria  Method 


The  Modal  Assurance  Criteria  (MAC)  and  the  method  of  sec¬ 
tion  2.4  utilize  theories  in  random  vibration  analysis.  For 
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TABLE  III 


EXPERIMENTAL  NATURAL  FREQUENCIES  FOR 
THE  SERIES  3  STABILATOR 


SINE  DWELL (cps) 
58 

101 

116 

124 

161 

198 

260 

279 

300 


MODAL  ASSURANCE 
CRITERIA(cps) 

56 

74 

101 

116 

124 

135 

146 

160 

176 

195 

258 

278 

280 

286 

300 


HP5451B  FOURIER 
ANALYSER(cps) 


102 

124 


160 

176 

198 

258 


286 
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discussions  of  the  theory  of  random  vibrations  the  reader  is 
referred  to  Ref  9,10.  The  primary  tool  for  both  investiga¬ 
tions  was  the  Hewlett  Packard  HP5451B  Fourier  Analyser  (Ref  11). 

The  author  was  introduced  to  the  MAC  test  through  the 
thesis  by  Larry  B.  Glenesk  (Ref  12)  as  an  accurate  means  of 
pinpointing  natural  frequencies  of  complex  structures.  A  dis¬ 
cussion  of  the  function  and  the  parameters  in  the  function  is 
presented  in  Ref  12:24-32.  The  MAC  function  is  a  measure  of 
the  coherence  between  responses  at  two  points  on  a  structure 
due  to  an  impulse  input  at  an  arbitrary  point  (striking  the 
structure).  The  MAC  function  is  defined  as: 


MAC 


JJL 


(co) 


S  (w)  S  (  u> ) 
yy  rr 


where  Sry(u)  represents  the  cross  power  spectral  density  of 
the  functions  * y f  and  'r'  in  the  frequency  (u>)  domain.  The 
functions  1 y *  and  'r'  would  represent  the  outputs  of  the  refer¬ 
ence  and  response  accelerometers  in  this  case.  The  bar  denotes 
an  ensemble  average  of  the  power  spectra  due  to  excitations  at 
any  number  of  different  points. 


2.3.2  Test  Setup 

Before  beginning  the  test,  the  stab  was  marked  with  an 
eleven  spanwise  by  eight  chordwise  grid  (Fig  2.2).  Table  IV 
is  a  list  of  the  grid  points  (CP's)  and  their  respective  HSS 
and  FS.  This  particular  grid  size  was  chosen  because  it  was 
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Figure  2.2,  Vibration  Test  Grid 


TABLE  IV 

SERIES  3  EXPERIMENTAL  GRID 


GRID 

POINT  HSS  FS 


1 

29.0 

544.5 

2 

29.0 

536.0 

3 

29.0 

527.5 

4 

29.0 

519.0 

5 

29.0 

515.3 

6 

29.0 

510.5 

7 

29.0 

502.0 

8 

29.0 

493.5 

9 

29.0 

485.0 

10 

33.3 

487.8 

11 

33.3 

495.8 

12 

33.3 

503.9 

13 

33.3 

512.0 

14 

33.3 

520.0 

15 

33.3 

528.1 

16 

33.3 

536.2 

17 

33.3 

544.2 

18 

39.1 

543.9 

19 

39.1 

536.4 

20 

39.1 

528.9 

21 

39.1 

521.4 

22 

39.1 

513.9 

23 

39.1 

506.4 

24 

39.1 

498.9 

25 

39.1 

491.4 

26 

44.8 

495.1 

27 

44.8 

502.0 

28 

44.8 

509.0 

29 

44.8 

515.9 

30 

44.8 

522.8 

31 

44.8 

529.7 

32 

44.8 

536.6 

33 

44.8 

543.5 

34 

50.5 

543.2 

35 

50.5 

536.8 

36 

50.5 

530.5 

37 

50.5 

524.2 

38 

50.5 

517.8 

39 

50.5 

511.5 

40 

50.5 

505.2 

41 

50.5 

498.8 

42 

56.3 

502.5 

43 

56.3 

508.3 

44 

56.3 

514.0 

45 

56.3 

519.8 

46 

56.3 

525.5 

GRID 

POINT  HSS  FS 


47 

56.3 

531.3 

48 

56.3 

537.0 

49 

56.3 

542.8 

50 

62.0 

542.4 

51 

62.0 

537.3 

52 

62.0 

532.1 

53 

62.0 

526.9 

54 

62.0 

521.7 

55 

62.0 

516.6 

56 

62.0 

511.4 

57 

62.0 

506.2 

58 

67.8 

509.9 

59 

67.8 

514.5 

60 

67.8 

519.1 

61 

67.8 

523.7 

62 

67.8 

528.3 

63 

67.8 

532.9 

64 

67.8 

537.5 

65 

67.8 

542.1 

66 

73.5 

541.7 

67 

73.5 

537.7 

68 

73.5 

533.7 

69 

73.5 

529.7 

70 

73.5 

525.6 

71 

73.5 

521.6 

72 

73.5 

517.6 

73 

73.5 

513.6 

74 

79.3 

517.3 

75 

79.3 

520.7 

76 

79.3 

524.2 

77 

79.3 

527.6 

78 

79.3 

531.0 

79 

79.3 

534.5 

80 

79.3 

537.9 

81 

79.3 

541.3 

82 

85.0 

541.0 

83 

85.0 

538.1 

84 

85.0 

535.3 

85 

85.0 

532.4 

86 

85.0 

529.5 

87 

85.0 

526.7 

88 

85.0 

523.8 

89 

85.0 

521.0 

90 

20.0 

515.3 

91 

10.0 

515.3 

92 

0.0 

515.3 
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fine  enough  to  identify  all  node  lines  from  0-300  cps. 

Two  accelerometers  were  placed  at  various  GP’s.  Their 
outputs  were  amplified  and  monitored  in  the  test  chamber  on 
a  dual  trace  oscilliscope .  Outputs  were  also  passed  to  a 
separate  room  containing  the  HP5451B  (HP,  Fig  2.3)  Fourier 
Analyser  where  they  were  filtered  to  below  500  cps  and 
digitized . 

2.3.3  Data  Aquisition 


One  accelerometer  remained  stationary  through  the  entirety 
of  each  testing  session  and  was  referred  to  as  the  reference 
accelerometer.  The  other  was  placed  at  arbitrary  points  on 
the  structure  and  was  referred  to  as  the  response  accelerometer 


Figure  2.3,  Hewlett  Packard  HP5451B  Fourier  Analyser 


To  discover  a  maximum  number  of  natural  frequencies  in  a 
test  session  the  reference  accelerometer  should  lie  at  a  point 
free  of  nodes  in  the  desired  frequency  range.  Following  sug¬ 
gestions  from  the  technical  staff  at  AFWAL/FIBG,  points  along 
the  outboard  trailing  edge  were  chosen. 

Table  V  is  a  list  of  reference  points  and  response  points 
tested.  Appendix  B,  program  1  presents  the  program  used  by 
the  HP  while  aquiring  data  and  processing  the  MAC  function. 

For  each  pair  of  points  the  responses  to  an  impulse  were 
recorded  and  averaged  in  an  attempt  to  reduce  the  affects  of 
random  noise.  At  the  end  of  each  ensemble  average  the  MAC 
was  calculated  and  displayed.  Judgment  of  its  accuracy  was 
based  on  the  number  of  spikes  that  peaked  between  zero  and 
unity.  A  large  number  of  spikes  appearing  say  between  10  and 
90  percent  of  unity  would  suggest  extensive  noise  corruption 
or  faulty  equipment.  Accurate  MAC  functions  were  produced  in 
hard  copy,  both  graphically  and  numerically,  for  each  pair  of 
points  on  a  Versatek  copier  (see  Fig  2.4  for  a  sample  graph). 
Some  pairs  required  several  attempts  to  acquire  accurate 
results . 

2.3.4  Results  of  the  MAC  Test 

Values  of  unity  on  the  graph  in  Fig  2.4  represent  natural 
frequencies.  Once  an  approximate  value  for  a  resonant  fre¬ 
quency  was  found,  the  data  was  viewed  numerically  so  that  a 
more  accurate  prediction  could  be  made.  When  successive 


TABLE  V 

REFERENCE  AND  RESPONSE  POINTS 
USED  IN  THE  MAC  TEST 


REFERENCE 

RESPONSE 

POINT 

POINT 

66 

15 

66 

31 

66 

47 

66 

63 

66 

79 

66 

11 

66 

27 

66 

43 

66 

59 

66 

75 

82 

12 

82 

19 

82 

20 

82 

57 

values  of  numbers  within  0.2  percent  of  unity  were  encountered , 
the  mean  frequency  was  recorded. 

Table  III  contains  the  natural  frequencies  found  using 
the  MAC  test.  Examination  of  Table  III  shows  that  the  MAC 
has  predicted  all  the  frequencies  found  in  the  Sine  Dwell  test 
plus  some  additional  ones. 

2 .4  Modal  Analysis  Using  the  HP5451B  Modal  Analysis  Software 
Package 


2.4.1  Overview 

To  make  a  complete  comparison  of  the  Series  3  stab  model 
required  a  visual  comparison  of  experimental  and  numerically 
predicted  mode  shapes  as  well  as  a  list  of  natural  frequencies. 
As  shown  by  Glenesk  (Ref  12),  experimental  mode  shapes  could 
be  generated  using  the  MAC  test  but  the  procedure  is  arduous 
and  time  consuming.  At  the  suggestion  of  Mr.  Richard  D. 
Talmadge,  AFWAL/FIBG,  it  was  decided  to  use  a  modal  analysis 
software  package  developed  for  use  on  the  HP  Fourier  Analysis 
System  (Ref  13). 

Experimental  modal  analysis  using  such  a  complex  system 
provides  enough  information  and  topics  for  discussion  to 
easily  render  a  thesis  by  itself.  Because  modal  analysis  is 
not  the  topic  of  this  thesis,  many  facts  concerning  the  test 
and  results  will  be  treated  in  a  'black  box'  fashion. 
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2.4.2  Experimental  Model  and  Test  Setup 


To  generate  a  mode  shape  for  the  stab,  the  HP  system 
required  information  on  model  geometry.  Reference  10,  sec¬ 
tion  3  discusses  the  procedure  and  information  required  to 
build  a  model.  Appendix  C  presents  this  information  for  the 
stab.  A  comparison  of  the  coordinates  listed  in  Appendix  C 
and  those  given  in  Table  IV  yields  the  respective  HSS  and  FS. 
Any  mention  of  grid  points  in  this  section  refers  specifically 
to  this  grid. 

In  contrast  to  the  MAC  test,  the  HP  system  requires  in¬ 
formation  about  the  impulse  to  the  structure  as  well  as  the 
response.  This  required  the  use  of  a  load  cell  hammer 
(Fig  2.5)  and  one  accelerometer.  With  this  exception,  the 
test  setup  was  identical  to  that  of  the  MAC  test. 


Figure  2.5,  Load  Cell  Hammer  Used 
in  Experimental  Modal  Analysis 
33 


2.4.3  Data  Acquisition 


Throughout  the  test  the  structure  was  excited  with  the 
load  cell  hammer  at  GP  7.  GP  7  was  chosen  because  of  its  loca¬ 
tion  along  the  root  rib.  Stiffness  of  the  root  rib,  as  opposed 
to  locations  on  the  structure  composed  only  of  skin  and  honey¬ 
comb,  prevented  structural  deformation  which  would  have  corrupted 
the  impulse.  In  addition,  GP  7  was  far  enough  from  the  elastic 
axis  to  excite  both  torsion  and  bending  modes. 

Five  samples  per  ensemble  per  GP  were  used  (the  structure 
was  struck  and  the  response  measured  five  times)  as  opposed  to 
fifteen  for  the  MAC  test.  If  the  test  went  perfectly,  the 
structure  would  have  been  struck  465  times.  This  illustrates 
the  reason  for  five  averages  as  opposed  to  the  fifteen  used 
for  the  MAC  test. 

One  impulse  and  response  ensemble  was  taken  for  each  GP. 

At  the  end  of  the  averaging  process  the  autocorrelation,  auto¬ 
power  spectral  density  (PSD),  cross  correlation,  cross  PSD,  and 
transfer  function  were  all  calculated.  At  this  time  a  decision 
was  made  as  to  the  accuracy  of  the  ensemble  by  evaluating  the 
above  information.  Accurate  transfer  functions  were  stored  by 
the  HP  for  use  in  data  reduction  after  the  completion  of  the 
test.  Appendix  B,  program  2  is  the  program  used  by  the  HP 
during  the  test. 

To  prevent  aliasing  of  the  input  information  the  response 
was  given  artificial  damping  by  multiplying  it  by  an  exponen¬ 
tial  decay  function  that  began  at  unity  and  became  zero  toware 
the  end  of  the  time  window.  The  impulse  was  handled  in  a 
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similar  manner  by  multiplying  it  by  a  truncated  sine  wave. 
This  prevented  errors  in  the  Fast  Fourier  Transforms. 


2.4.4  Data  Reduction  and  Results 


To  generate  a  number  of  mode  shapes  the  HP  system  first 
required  a  list  of  frequencies  of  interest.  These  could  be 
entered  manually  or  picked  from  a  representative  transfer  func¬ 
tion  using  a  display  screen  and  cursor.  The  latter  was  used 
based  on  the  fact  that  not  all  natural  frequencies  found  using 
the  earlier  methods  were  prominent  in  the  majority  of  transfer 
functions . 

Once  a  list  of  frequencies  had  been  picked,  a  method  of 
curve  fitting  was  chosen.  The  methods  available  are  presented 
in  Ref  13,  section  4.  Because  damping  information  was  not 
desired,  the  Quadrature  method  was  used  to  find  a  relative 
amplitude  for  each  point  at  each  of  the  chosen  frequencies. 
Results  of  the  Quadrature  fit  for  a  particular  frequency  could 
be  previewed  in  tabular  form  along  with  an  animated  display  of 
the  mode  shape.  Figures  2.6  and  2.7  are  representative  of  the 
modes  extracted  using  this  technique.  Appendix  D  contains 
pictures  of  all  the  modes  found  between  0-400  cps.  Some  care 
should  be  exercised  in  viewing  the  modes  above  260  cps  because 
the  grid  may  not  have  been  dense  enough  to  show  all  the  node 
lines  present.  The  superposition  of  a  number  of  displacement 
amplitudes  through  a  complete  cycle  highlights  the  node  lines. 
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Experimental  Mode  Shape,  100  cps 


Lgure  2.7,  Experimental  Mode  Shape,  124  cps 


2.4.5  Comments 


Some  question  may  arise  as  to  why  the  procedure  discussed 
in  this  section  did  not  yield  mode  shapes  for  all  the  natural 
frequencies  predicted.  The  answer  lies  in  the  inability  of 
the  curve  fitting  technique  to  differentiate  between  modes 
lying  within  ±  3  cps  of  each  other.  Another  reason  lies  in 
the  fact  that  transfer  functions  in  the  frequency  domain  are 
representative  of  the  amount  of  energy  in  a  mode  of  vibration, 
i.e.  low  energy  modes  have  small  peaks  and  are  hard  to  fit 
numerically. 

The  first  reason  discussed  above  is  the  most  significant 
reason  for  not  extracting  a  mode  shape  at  58  cps.  Because 
60  cps  line  noise  was  prominent  in  the  lab  it  could  not  be 
averaged  out  and  consistently  foiled  attempts  to  fit  the 
58  cps  mode.  The  74  cps  mode  predicted  by  the  MAC  test  could 
not  be  fit  because  its  magnitude  was  consistently  smaller  than 
the  60  cps  noise  peak.  Attempts  were  made  to  increase  the 
magnitude  by  introducing  additional  energy  in  the  impulse  and 
low  pass  banding  the  response  to  below  200  cps.  A  new  rubber 
tipped  hammer  was  designed  (Fig  2.8)  with  an  additional  lumped 
mass  above  it  to  accomplish  this.  The  attempt  failed  because 
the  automatic  scaling  factor  in  the  HP  was  set  each  time  by 
the  60  cps  noise  peak.  This  made  the  74  cps  mode  consistently 
appear  diminutive. 
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Figure  2.8,  Modified  Load  Cell  Hammer 

2 . 5  Numerical  Eigenvalue  Analysis  of  the  Series  2  and 
and  Series _ 3  Finite  Element  Models 


2.5.1  Analysis  Method 


With  the  experimental  results  obtained  from  the  methods 
discussed  in  sections  2. 2-2. 4,  an  accurate  appraisal  of  the 
stiffness  distribution  of  the  Series  3  model  can  be  made. 
NASTRAN  Rigid  Format  3  was  vised  to  numerically  predict 
natural  frequencies  and  mode  shapes  for  the  Series  3  stab  in 
a  Free-Free  condition.  The  inverse  power  method  was  used  with 
normalization  accomplished  on  the  largest  value  in  the  analysis 
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set  (Ref  7).  The  NASTRAN  1 SUPORT’  card  was  used  to  suppress 
the  rigid  body  modes. 

To  further  investigate  the  statement  made  in  Ref  5  that 
the  Series  2  and  Series  3  stabs  have  equivalent  stiffness,  the 
natural  frequencies  of  the  Series  2  stab  were  also  predicted. 

2.5.2  Results  and  Conclusions  for  the  Eigenvalue  Analysis 

Table  VI  contains  the  list  of  natural  frequencies  pre¬ 
dicted  numerically  for  the  Series  2  and  Series  3  stabs  from 
0-300  cps.  A  comparison  shows  them  to  be  virtually  identical. 
These  results  suggest  dynamic  similarity  while  statically  the 
structures  react  differently  as  was  shown  in  Chapter  1.  In 
this  respect,  the  statement  that  the  two  possess  equivalent 
stiffness  is  true. 

Mode  shapes  for  the  Series  3  stab  were  generated  on  a 
Tektronix  4014  terminal  using  the  interactive  program 
GCSNAST  (Ref  11).  Figures  2.9  and  2.10  are  representative 
mode  shapes.  Appendix  D  contains  mode  shapes  numerically 
predicted  for  the  Series  3  stab  up  to  400  cps.  Some  care 
should  be  exercised  when  viewing  mode  shapes  above  300  cps 
because  the  grid  may  not  be  dense  enough  to  show  all  the  node 
1 ines . 
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TABLE  VI 

CALCULATED  NATURAL  FREQUENCIES 


SERIES  2 

NATURAL  FREQUENCY(cps) 

SERIES  3 

NATURAL  FREQUENCY (cps) 

46.45 

46.65 

78.66 

77.70 

88.31 

88.04 

115.39 

112.97 

122.56 

121.56 

126.76 

124.76 

147.35 

144.85 

168.64 

167.09 

189.54 

188.73 

199.91 

200.05 

213.82 

214.32 

224.25 

224.78 

252.81 

252.46 

258.83 

259.82 

280.64 

279.62 

310.94 

309.24 
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2.6  Conclusions  From  Modal  Analysis 


Figures  2.6,  2.9  and  2.7,  2.10  are  experimental  and 
theoretical  mode  shapes  that  appear  to  match  exactly.  Others 
that  show  a  marked  resemblance  are  presented  in  Appendix  D 
and  their  relationship  is  shown  in  Table  VII. 


TABLE  VII 

EXPERIMENTAL  AND  NUMERICALLY  PREDICTED 
MODE  SHAPES  THAT  MATCH  IN  APPENDIX  D 


EXPERIMENTAL 

Fig 

Freq  (cps) 

D1 

100 

D2 

124 

D3 

160 

D4 

176 

D5 

198 

D6 

258 

D7 

286 

D8 

352 

NUMERICAL 

Fig 

Freq  (cps) 

D12 

77.7 

D13 

88.0 

D14 

113.0 

D15 

121.6 

D18 

167.1 

D21 

214.3 

D22 

224.8 

D24 

259.8 

While  NASTRAN  has  predicted  the  shape  accurately,  the  fre¬ 
quencies  predicted  are  20-40  cps  too  low  in  the  majority  of 
the  cases.  This  suggests  errors  in  the  stiffness  of  the  model 
due  to  either  the  mass  distribution  or  material  stiffness 
(elastic  modulus). 

As  pointed  out  in  Chapter  1,  the  Series  3  model,  was  found 
to  be  lacking  in  stiffness  when  subject  to  static  loads.  When 
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Figure  2 .10 , Numerically  Predicted  Mode  Shape,  88  cps 


the  relationship  between  frequency  and  elastic  modulus  for  a 
simple  beam  is  considered,  i.e.  frequency  is  proportional  to 
the  square  root  of  Young's  modulus,  assuming  no  error  in 
moment  of  inertia,  it  is  logical  from  the  static  results  that 
the  predicted  frequencies  would  be  too  low.  This  relationship 
can  be  used  to  make  a  crude  approximation  to  the  change  in 
overall  stiffness  required  to  'tune'  the  model  into  performing 
correctly. 

Because  finite  element  models  contain  finite  numbers  of 
nodes,  it  is  expected  that  the  error  involved  in  predicting 
natural  frequencies  gets  worse  as  the  frequency  increases. 

This  is  exhibited  by  the  fact  that  the  differences  in  frequency 
of  the  modes  predicted  correctly  is  not  a  constant.  For  this 
reason,  the  first  experimental  and  numerically  predicted  mode 
shapes  that  matched  were  used  to  predict  a  correction 
( see  Figs  2.6,  2.9). 


E  =  E  ,  ,X 
new  old 


100 

88 


=  1.64  E 


old 


This  suggests  an  overall  increase  in  stiffness  of  64  percent. 

Believing  this  to  be  somewhat  extreme  it  was  decided  to  exper¬ 
iment  with  a  value  of  30  percent. 

Chapter  1,  section  1.3  provides  a  list  of  the  elements  of 
the  Series  3  stab  affected  by  an  increase  in  shear  and  elastic 
moduli  in  an  attempt  to  correct  the  model.  An  eigenvalue 
analysis  of  the  stiffened  Series  3  model  was  done  but  time 
constraints  and  limited  file  space  on  the  AFIT  computer  system 
prevented  generation  of  new  mode  shapes.  However,  a  comparison 
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of  the  displacements  of  GP's  in  a  normalized  mode  shape  listed 
in  the  NASTRAN  output  showed  a  shift  of  10  cps  upward  for  each 
of  the  two  modes  in  figures  2.9  and  2.10.  These  were  the  only 
modes  identifiable  in  this  manner. 

Reference  4,  section  III  provides  a  comparison  of  the 
mass  distribution  of  the  model  with  information  from  Northrop 
report  NAI-58-11  (Ref  15).  The  correlation  is  good  with  the 
exception  of  50  percent  additional  mass  in  the  root  section. 

A  comparison  of  experimental  mode  shapes  immediately  after 
124  cps  and  theoretical  mode  shapes  after  88  cps  shows  similar 
node  patterns  except  for  an  extra  node  appearing  at  the  inboard 
trailing  edge  of  the  theoretical  model.  These  findings  suggest 
a  problem  in  either  the  mass  or  the  stiffness  of  the  elements 
in  this  area. 
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CHAPTER  3 


VERIFICATION  OF  A  STEADY  AERODYNAMIC 
MODEL  OF  THE  T-38  STABILATOR 


3.1  Introduction 


Building  a  stabilator  model  for  use  in  flutter  analysis 
requires  consideration  of  the  aerodynamic  loads.  Ultimately, 
verification  of  an  unsteady  aerodynamic  model  is  necessary. 

At  present  the  procedure  for  verifying  an  unsteady  model  has 
not  been  established  and  no  experimental  results  exist.  To 
add  some  confidence  in  the  model,  the  steady  aerodynamics  will 
be  considered. 

3 . 2  Analytical  Models 


Two  models  of  the  stab  were  constructed  for  analysis  on 
USSAERO  (Ref  16)  and  NASTRAN.  Appendix  E  contains  listings 
of  the  input  required  for  each  model.  Centerline  symmetry  was 
assumed  to  simulate  attachment  of  the  root  to  a  wall,  see 
Milne-Thomson  (Ref  17).  Neglecting  viscous  affects,  this  is 
an  accurate  representation  of  an  airfoil  in  a  large  diameter 
wind  tunnel. 

Both  the  USSAERO  and  NASTRAN  models  were  designed  to 


correlate  with  experimental  results  obtained  at  the 
University  of  Texas  at  Austin  (Ref  18).  Figure  3.1  (Ref  18) 
illustrates  the  geometry  of  the  experimental  model. 


3.2.1  USSAERO  Model 


For  an  extensive  discussion  of  the  aerodynamic  theory 
contained  in  USSAERO  the  reader  is  referred  directly  to  Ref  16. 
In  summary,  USSAERO  lifting  panels  are  modeled  using  a  smeared 
vorticity  distribution  varying  linearly  in  the  streamwise 
direction.  Singularity  strengths  are  related  to  normal  velo¬ 
city  by  satisfying  the  no  flow  condition  at  the  geometric  cen¬ 
ter  of  each  panel.  This  provides  a  system  of  equations  that 
yield  the  strength  of  the  singularity  distribution  over  each 
panel.  The  velocity  and  pressure  coefficients  are  then 
obtained  from  the  singularity  strengths.  The  aerodynamic 
force  on  a  panel  is  calculated  by  multiplying  the  pressure 
coefficient  by  the  panel  area. 

Two  options  are  available  for  modeling  thickness  affects. 
Either  the  upper  and  lower  surfaces  of  the  airfoil  can  be 
paneled  or  just  the  mean  camber  line.  If  the  mean  camber  line 
is  paneled,  sources  and  sinks  are  distributed  along  it  to  sim¬ 
ulate  thickness.  This  method  reduces  the  number  of  control 
points  and  the  calculation  time  considerably  and  was  chosen 
for  the  analysis.  No  consideration  for  anhedral  has  been  made 
because  of  the  small  angles  of  attack  tested. 
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Airfoil  section  -  65A004 

Planform  area  -  1.626  s 


Ref  18 


3.2.2  NASTRAN  Model 


Work  on  the  initial  development  of  the  aerodynamic  theory 
in  NASTRAN  is  presented  in  Ref  19,20.  Several  options  are 
available  depending  on  the  flight  regime  of  interest.  The 
Doublet  Lattice  Method  is  used  for  the  analysis  in  this  chap¬ 
ter  and  Chapter  4  and  is  summarized  in  the  following  paragraphs  . 
It  applies  only  to  subsonic  flight. 

NASTRAN  unsteady  lifting  panels  are  modeled  using  both 
a  horseshoe  vortex  and  acceleration  line  doublet.  The  horse¬ 
shoe  vortex  passes  through  the  quarter  chord  of  each  panel 
with  trailing  vorticies  entering  and  leaving  along  the  panel 
edges.  The  acceleration  line  doublet  lies  directly  on  the 
bound  part  of  the  horseshoe  vortex.  A  system  of  equations  is 
generated  by  relating  the  singularity  strengths  to  the  normal 
velocity  at  the  3/4  chord,  spanwise  center  of  each  panel  and 
staisfving  the  no  flow  condition.  Once  the  singularity 
strength  is  known,  the  Kutta-Joukowski  theorem  is  used  to  find 
the  magnitude  of  the  lift  at  the  1/4  chord  of  each  panel.  No 
account  is  made  of  thickness  or  anhedral  effects. 

The  NASTRAN  model  is  made  up  of  a  structural  element, 
'CQUAD1',  and  ninety-nine  aero-boxes,  'CAFR01'.  For  the  pro¬ 
cedure  used  to  generate  and  output  the  pressure  coefficient 
on  each  box,  the  reader  is  referred  to  the  thesis  by  Lance  E. 
Chrisenger,  AFIT/GAE/AA/80D-2  (Ref  21).  It  involves  forcing 
the  model  into  a  rigid  pitch  at  the  desired  angle  of  attack 


50 


and  approximately  zero  reduced  frequency.  Each  element  in  the 


'aeroforce'  vector  is  divided  by  its  respective  panel  area 
(constant  pressure  panel)  and  output  to  yield  values  of  the 
pressure  difference  at  the  1/4  chord. 


3 . 3  Method  of  Analysis  and  Results 


Experimental  information  from  the  University  of  Austin 
allowed  comparison  of  the  chordwise  pressure  distribution  on 
the  upper  surface  of  the  airfoil  at  two  angles  of  attack  at 
Mach  0.19.  Because  USSAERO  is  a  program  specifically  designed 
for  steady  aerodynamic  analysis,  pressure  distributions  are 
available  as  output  for  both  the  upper  and  lower  surfaces  of 
the  airfoil.  NASTRAN  however,  computes  resultant  forces  for 
use  in  flutter  analysis  and  only  the  pressure  difference  is 
available . 

Verification  of  the  NASTRAN  steady  aerodynamic  model  was 
accomplished  by  first  comparing  the  USSAERO  solution  to  the 
experimental  data  for  the  upper  surface  (Fig  3.2).  After 
verifying  the  USSAERO  model,  a  comparison  of  the  pressure  dif¬ 
ferences  predicted  by  USSAERO  and  NASTRAN  was  made  (Fig  3.3). 
Appendix  F  contains  comparisons  of  theoretical  and  experimen¬ 
tal  results  at  33.5,  66.9  and  91.7  percent  semi  span  for  two 
angles  of  attack  at  Mach  0.19. 
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66.9  PERCENT  SEMISPflN 


0 -USSAERO 


□  -  EXPERIMENTAL 
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Figure  3.2,  Comparison  of  Experimental  and 

Numerically  Predicted  Pressure  Distributions 
On  the  Stab  Upper  Surface,  a  -  2.3  M  =  0.19 
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66.9  PERCENT  SEMISPAN 


Figure  3.3,  Comparison  of  Numerical Iv  Predicted 
Pressure  Differences  on  the  Stab 


\k 


I 

j 

3 . 4  Conclusions 

Relatively  good  correlation  has  been  achieved  between 
USSAERO  and  experimental  data.  Small  differences  lie  along 
the  first  fifteen  percent  and  last  thirty  percent  of  the  chord. 

USSAERO  solutions  tended  toward  infinity  at  the  leading  edge 
as  is  inherent  in  potential  aerodynamic  techniques.  In  each 
case  the  numerical  solution  tends  to  intersect  the  experimen¬ 
tal  solution  but  not  at  any  specific  point.  This  suggests  an 
error  in  the  moment  generated  over  the  airfoil  even  though  an 
accurate  value  of  section  lift  may  be  predicted.  One  solu¬ 
tion  mav  be  to  concentrate  panels  toward  the  leading  edge  of 
the  airfoil  to  closer  aporoximate  a  curved  pressure  rise  in 
that  area. 

Comparison  of  pressure  differences  across  the  airfoil 
using  NASTRAN  and  USSAERO  show  extremely  close  results. 

USSAERO  has  predicted  consistently  higher  pressure  differences 
than  NASTRAN  but  the  differences  are  small  and  not  cause  for 
concern.  This  was  expected  due  to  the  lack  of  thickness 
affects  in  the  NASTRAN  model.  These  results  show  that 
NASTRAN  is  modeling  the  aerodynamics  of  the  airfoil  correctly 
for  the  steadv  case. 
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CHAPTER  4 


FLUTTER  ANALYSIS  OF  THE  T-38 
STABILATOR  USING  NASTRAN 


4.1  Introduction 


Accurate  modes  of  vibration  and  natural  frequencies  are 
of  extreme  importance  in  determining  flutter  speed.  Under 
different  flight  conditions  an  airfoil  may  flutter  in  any  mode 
but  most  often  flutter  is  encountered  in  one  of  the  first 
modes,  i.e.  torsion,  bending,  or  a  combination  of  the  two.  In 
Ref  4,  "three  and  possibly  four  modes"  were  found  that  corres¬ 
pond  closely  to  those  of  the  stab  in  a  flexible  root 
(1  hvdraulic  actuator  system  on)  boundary  condition.  Resting 
on  these  results,  a  method  of  more  accurately  representing  the 
aerodynamics  of  the  stab  in  an  in-flight  condition  was  consid¬ 
ered. 


4.2  Analvtical  Model 


The  NASTRAN  flutter  model  is  composed  of  the  right  and 
left  stabilators  and  a  mock  fuselage  but  no  wing.  In  Ref  15 
the  wing  is  said  to  be  stable  up  to  115  percent  of  the  flight 
envelope.  Perturbations  to  the  wing  may  cause  it  to  oscillate 
but  within  the  flight  envelope  the  motion  damps  out.  At  some 
frequency  of  oscillation  it  is  probable  that  flutter  in  the 
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stab  may  be  excited  by  the  wing  wake.  Although  this  may  be  a 
problem,  no  attempt  was  made  to  include  effects  from  the  wing 
because  of  the  complex  structural  model  required. 

The  stab  is  divided  into  six  spanwise  and  four  chordwise 
panels.  This  rather  sparse  panelling  scheme  was  chosen  follow¬ 
ing  discussions  with  engineers  at  AFWAL/FDL  who  had  achieved 
good  results  with  similar  models.  A  change  in  the  geometry  at 
the  root  of  the  stabilizer  was  made  to  straighten  it.  This 
was  required  in  order  to  panel  the  stab  and  insure  that  panel 
side  edges  be  parallel  to  the  free  stream.  It  involved  a 
minor  shift  of  grid  points  136-144  from  their  respective  HSS 
to  HSS  27,567. 

NASTRAN  'slender  bodies'  and  'interference  bodies'  are 
used  to  model  the  fuselage  from  the  wing  apex  to  the  stabila- 
tor  inboard  trailing  edge  (Ref  7,20).  Slender  bodies  are 
composed  of  doublets  placed  in  the  flow  field  to  simulate 
fulselage  thickness.  A  circular  or  elliptical  fuselage  cross 
section  can  be  generated  by  placing  singularities  at  the  cen¬ 
terline  or  equidistant  from  the  centerline  in  the  aircraft 
X-Y  plane.  Interference  bodies  are  circular. 

Elliptical  slender  bodies  were  chosen  with  major  and 
minor  axes  corresponding  to  the  width  and  height  of  the  cross 
section  at  the  stab  attachment  point,  FS  516.  Eleven  slender 
bodies  are  used  to  satisfy  a  requirement  made  in  Ref  20, 

Part  II  that  the  ratio  of  major  axis  to  length  be  equivalent 
or  greater  than  two.  This  assures  a  number  of  doublets  suf¬ 
ficient  to  simulate  at  least  a  flat,  continuous  surface.  The 
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fuselage  cross  section  area  is  reduced  linearly  to  zero  from 
the  wing  1/4  chord  to  the  wing  apex  to  simulate,  in  affect, 
a  conical  nose. 

Interference  bodies  use  the  method  of  images  to  satisfy 
the  no  flow  condition  through  the  fuselage.  Reference  20 
suggests  the  possibility  of  creating  interference  bodies  of 
noncircular  cross  section  but  only  circular  cross  sections  are 
available  in  NASTRAN.  The  method  of  images  involves  placing 
singularities  within  the  radius  of  the  interference  body  to 
negate  the  normal  wash  on  its  surface.  Three  interference 
bodies  are  used  with  division  points  at  the  wing  apex,  stab 
inboard  leading  edge,  midroot  chord,  and  trailing  edge.  A 
constant  radius  of  27.567  inches  (the  half  width  of  the  fuse¬ 
lage  at  the  stab  attachment  point,  FS  516)  was  used. 

4. 3  Method  of  Analysis 

Flutter  analysis  using  NASTRAN  is  accomplished  with 
Rigid  Format  10.  The  method  of  analysis  is  designated  on  the 
NASTRAN  'FLUTTER'  card  as  the  'K'  method  (Ref  22,  section  17). 
The  ’ K '  method  involves  looping  through  values  of  density 
ratio,  reduced  frequency,  and  Mach  number.  Values  of  flutter 
velocitv,  damping  ratio,  and  flutter  frequency  are  calculated 
at  the  end  of  each  iteration  and  presented  in  the  output. 

Prior  to  considering  the  flutter  problem,  the  aerodynamic 
modal  force  vector  Qhh  must  be  generated.  Elements  in  the 


vector  are  functions  of  reduced  frequency  (k)  and  Mach  number 
(M)  that  are  specified  on  the  'MKEAR01'  card.  The  vector  is 
segmented  in  blocks  which  represent  specific  M,k  pairs.  For¬ 
mulation  of  Qhh  is  an  expensive  process  and  a  minimum  number 
of  Mach  numbers  and  reduced  frequencies  should  be  chosen. 
Obviously,  some  prior  knowledge  of  the  flutter  parameters  is 
desirable  when  making  the  selection. 

Values  of  density  ratio  refer  to  a  reference  density 
specified  on  the  'AERO'  card.  The  corresponding  speed  of 
sound  is  also  required  on  this  card. 

At  sea  level: 

p  =  1.147  x  10- 7  lbp  sec2  /  in^ 

a  =  1.34  x  10^  in/sec 

The  units  on  density  should  be  noted.  These  are  consistent 
with  the  use  of  pounds  force/square  inch  for  elastic  and  shear 
moduli . 

Looping  is  controlled  by  the  ’FLUTTER'  card.  The 
'FLUTTER'  card  contains  labels  for  'FLFACT'  cards  that  specify 
looping  parameters.  Any  of  the  values  specified  on  the 
'MKAER01 '  card  plus  additional  intermediate  values  can  be  used. 
If  the  aerodynamic  modal  force  for  a  specific  M,  k  pair  is  re¬ 
quested,  and  it  has  not  been  generated  via  the  'MKAER01 '  card, 
two  methods  of  interpolation  are  available  to  determine  it. 
These  are  specified  on  the  'FLUTTER'  card  as  'L'  or  'S', 
linear  and  surface  respectively.  If  more  than  one  Mach  num¬ 
ber  is  used  to  generate  the  aerodynamic  modal  force  vector 
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(it  is  customary  to  choose  several  values  of  k) ,  surface 
interpolation  is  used.  Otherwise,  linear  interpolation  must 
be  used. 

The  connection  between  the  aerodynamic  and  structural  de¬ 
grees  of  freedom  is  accomplished  using  splines.  For  the 
2-dimensional  model,  surface  splines  are  used.  The  ’ SPLINE1' 
card  is  used  to  designate  points  in  the  grid  to  be  splined. 

To  reduce  computation  time,  fifteen  points  on  the  structure 
were  splined.  No  guidance  is  available  on  the  actual  number 
of  points  necessary  and  it  was  thought  that  this  number  would 
be  sufficient.  The  grid  points  are  illustrated  schematically. 


Use  of  more  GP’s  has  been  found  to  lead  to  excessive  computa¬ 
tion  time  and  core  requirements  making  turn  around  time  imprac¬ 
tical.  See  Fig  1.4  for  the  precise  locations. 

Output  is  presented  in  tabular  form  with  the  option  of 
generating  the  classical  V-g  diagram.  Appendix  G  contains  a 
sample  in  tabular  form.  Sign  changes  in  the  damping  ratio  (g) 
represent  points  of  encipient  flutter.  Linear  interpolation 
is  then  used  to  pinpoint  values  of  flutter  speed  and  cyclic 
and  reduced  frequency.  If  no  flutter  condition  is  encount¬ 
ered,  linear  extrapolation  can  be  used  to  provide  an  educated 
guess  at  new  values  of  M  and  k  to  begin  the  process  over. 
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4.4  Results  and  Conclusions 


Results  published  in  Ref  2  Cor  Clutter  velocity  oC  the 
stab  at  sea  level  were  used  as  initial  parameters  for  the 
Clutter  analysis  on  NASTRAN.  ReCerence  2  uses  Strip  Theory 
with  a  rigid  fuselage.  It  is  not  stated  in  ReC  2  as  to 
whether  wing  effects  were  considered.  Results  are  shown  in 
Table  VIII. 


TABLE  VIII 

COMPARISON  OF  NUMERICALLY  PREDICTED 
FLUTTER  SPEEDS 


STRIP  THEORY 
(Ref  2) 

NASTRAN 

Flutter  Speed  (KEAS) 

481.80 

494 

Frequency  (cps) 

30.51 

23.8 

Mach  Number 

0.73 

0.75 

The  Clutter  speeds  are  very  close  although  the  frequency  is 
slightly  off. 

This  indicates  good  correlation  between  the  two  methods 
with  NASTRAN  predicting  a  higher  flutter  speed.  These  results 
were  achieved  with  no  corrections  in  stiffness  of  the  model. 

It  is  expected  that  much  less  conservative  values  would  be 
calculated  using  a  corrected  model. 
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SUMMARY  OF  RESULTS  AND  CONCLUSIONS 


In  Chapter  1,  static  analysis  of  the  Series  3  finite 
element  model  revealed  considerable  lack  of  stiffness  in  chord- 
wise  bending.  This  conclusion  was  arrived  at  through  a  com¬ 
parison  of  the  characteristics  of  the  Series  2  finite  element 
model  with  experimental  results.  The  Series  2  and  Series  3 
models  are  identical  with  the  exception  of  three  intermediate 
ribs  and  an  auxilliary  spar  included  in  the  former. 

The  Series  2  model  was  found  to  lack  approximately  30 
percent  stiffness  in  chordwise  bending.  This  may  be  attribu¬ 
ted  to  modeling  the  structure  with  flat  plates  as  opposed  to 
torsion  cells  that  can  support  shear  flow.  Under  identical 
loads  and  boundary  conditions,  deflections  of  the  Series  3 
model  are  consistently  greater  than  those  of  the  Series  2 
model.  This  suggests  an  appreciable  difference  in  the  stiff¬ 
ness  of  the  models  and  refutes  the  statement  made  in  Ref  5 
that  they  are  equivalent.  Because  the  structures  are  virtu¬ 
ally  identical,  and  because  there  are  no  experimental  static 
test  results  for  the  Series  3  stab,  it  is  assumed  that  the 
Series  3  model  is  lacking  in  stiffness  by  approximately  30 
percent  also. 

Results  of  Chapter  2  support  the  hypothesis  that  the 
Series  3  model  is  lacking  in  stiffness.  A  numerical  eigen¬ 
value  analysis  of  the  model  yielded  natural  frequencies  con¬ 
sistently  lower  than  those  suggested  from  experiment . 
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Investigation  of  the  higher  order  coupled  modes  revealed  a 
problem  area  at  the  root,  trailing  edge.  It  was  shown  in 
Ref  4  that  the  mass  distribution  of  the  model  compared  well 
with  the  actual  stab  except  for  a  50  percent  deviation  at  the 
root.  The  error  is  attributed  to  this.  Experimental  mode 
shapes  at  100  cps  and  124  cps  were  matched  with  numerically 
generated  shapes  at  78  cps  and  88  cps.  Other  shapes  show 
marked  resemblance  all  the  way  to  350  cps  and  are  identified 
in  Table  VII.  The  deviation  in  shape  is  also  attributed  to 
the  error  discussed  above. 

Investigation  of  the  pressure  distribution  over  the  stab 
in  a  steady  condition  showed  good  results.  Integration  of  the 
chordwise  pressure  distribution  provides  accurate  section  lift 
coefficients.  However,  some  small  deviation  in  section  moment 
exists.  This  can  be  attributed  to  boundary  layer  affects  in 
the  experimental  results. 

Errors  in  model  stiffness  identified  in  Chapters  1  and  2 
suggest  predictions  of  flutter  speed  using  this  model  would  be 
conservative.  Results  of  Chapter  4  at  sea  level  show  a 
flutter  speed  equivalent  to  that  predicted  in  Ref  2  where  the 
method  of  analysis  was  Strip  Theory.  It  is  expected  that  a 
corrected  model  would  predict  flutter  speeds  much  higher  than 
those  found  using  Strip  Theory. 
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Experimental  measurement  of  the  displacements  of  the 
Series  3  stab  under  static  loads  is  desired.  With  these 
results  a  more  accurate  judgment  of  rhe  characteristics  of 
the  Series  3  model  and  corrections  can  be  made. 

No  investigation  of  the  difference  between  the  published 
section  mass  of  the  root  (Ref  15)  and  the  calculated  mass  of 
the  model  at  the  root  (Ref  4)  has  been  attempted.  Since  the 
two  differ  by  an  appreciable  amount,  and  a  problem  has  been 
identified  in  the  immediate  area,  some  investigation  is 
necessary. 

In  addition  it  is  desirable  to  repeat  the  experimental 
modal  analysis  discussed  in  Chapter  2,  section  2.4  for  bound¬ 
ary  conditions  corresponding  to  the  stab  attached  to  the  air¬ 
craft.  This  would  provide  first  hand  verification  of  the 
model  in  a  flutter  condition. 

Once  the  structural  aspects  of  the  model  are  verified, 
the  torsional  spring  stiffness  of  the  hydraulic  actuator  sy¬ 
stem,  noted  in  Ref  14  as  "critical”  to  the  flutter  speed, 
should  be  investigated.  Then  an  accurate  assessment  of  the 
model's  capabilities  in  flutter  speed  prediction  can  be  made. 
At  this  point,  a  procedure  for  simulating  damage  and  repair 
can  be  investigated. 
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APPENDIX  A 

COMPARISON  OF  DEFLECTIONS  OF  THE  SERIES  2 
AND  STIFFENED  SERIES  2  MODELS,  UNDER 
IDENTICAL  STATIC  LOADS, 

WITH  EXPERIMENTAL  RESULTS 
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Blacked  in  symbols  are  Series  2  results. 

Open  symbols  are  N0R-60-6  experimental  results. 
Dotted  lines  are  Stiffened  Series  2  results. 


LOAD (lb) 

HSS 

%  CHORD 

FIGURE 

400 

82.0 

Leading  edge 

A1 

400 

70.75 

20.0 

A2 

400 

70.75 

75.0 

A3 

1200 

42.25 

20.0 

A4 

1200 

42.25 

75.0 

A  5 

^  Load  applied  here. 
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Leading  Edge 
20%  Chord 
52.71  Chord 
75%  Chord 

Trailing  Edge 


Percent  Chord  at  HSS  30.00 


•  •  « 

I  •  .  • 

(sotpui)  uoi^osxjaa 
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Figure  A2 


Percent  Chord  at  HSS  30.0 


Percent 


APPENDIX  B 

PROGRAMS  USED  BY  THE  HP5451B 
FOURIER  ANALYSER 


73 


r-®T-4P)© 

*H(U 

coco  oo 
inin  in 


>>e>>-o  . 

njinoO’H^rt^ 

ojcucumnn 


© 

© 

in 


fu©cu 


<s> 

o 

in 


©  ©  © 
©  © 
t'"  -r-l 


COCO 


*-i  ©  ©  ru  ui  ru  © 


©  ©  ©  ^-# 
©® 
r- 1~> 


© 


^rj 
CU  H 
CO© 

inn 


(UHTUirt 


© 


»-)©©*-<©  'r-f  © 


H“r(6n 
■m  cu 
co  co 
in  in 


tun^-ru'Kj.'Sftyira 

*"*  ©  ©  O  r-t  ClJ  — 4 

CO  00  CO  OO  CO  CO  CO  00  CO 

mm  mm in  in  in m  in 


►J.. 

*-•>-> 


*0  . — )  . — i  ,  1  dT  W  r) 

o  o  o  o  m  ocnix 


JJO. 

OOtPtt 


I  VJ 

OQ>3>-JCI- J 


>>>>>->->• 


T  CO  CU  LD  CO  *-(  T  t'-  CO  CU  *r  >  -m  in  C:  T  O' 

t*t4  «-<(u£tx(Ucunn  co  in 


to-  ©  m  w  -i  co 

COdir-C-C-lJCO 


to  ©  ro  co  n  in 
£0  Or  cc-  O-  ©  © 


0>ino> 


74 


min  to 


OlvtU  OCx2C0'^'S>O)'r-lt^— t  -rHfU 

oo  ro^runj 

00 
in 


into 

QX>  o*i  =  si>-!r«>:^3caci-3  * 

in 

•e-  r -  cb  -<  n-  in  nj  ’S)  -~'T  r-  <9  no  to  o>  h  t  r- 

_i_.r\jiDnjnjr»rTr  —  Ti/unininiaiotfl 

'  —4  "«—*  W  — H  «— (HHH- 


•Si 


tu  ©  ni 


co  »-<  n-  «-•  ^  rj  o  <s>  ■-*  © 

m  co  «•«  m— *  ru 

© 


vnv\  ovv>  c>  tr>  jjj-j  c 

JECCCXr.ECXiJOOUOJ>E«X1tXOuJLOOIf 

t«  ■*•£•'■  o  ©  ro  f^-srotoer'WLooo-*  '4-30®  n  os  ct  ru  *rr-  ®  ■4-co  ru 
^•^•HnjAjfonnPin'n-j  u*u  u<u.'tOv£>ui('  i- (•-  wfufopi 


37  /L  100 
109  i  1 
195  V  301S 
109  CH  0 


APPENDIX  C 
THE  HP  MODEL 
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EXPERIMENTAL  AND  NUMERICALLY 
PREDICTED  MODE  SHAPES 
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Figure  D1 ,  100  cps 


Figure  D3,  160  cps 


AIR  FORCE  INST  OF  TECH  MR  1 GHT-P A  TTERSON  AFB  OH  SCHOO— ETC  F/6  1/3 
INVESTIGATION  OF  AN  IMPROVED  FLUTTER  SPEED  PREDICTION  TECHNIQUE— ETC (U) 
DEC  80  R  K  THOMSON 

AFIT/GAE/AA/80D-21  NL 


Figure  D5,  198  cps 


Figure  D6,  258  cps 


Figure  D7,  286  cps 


Figure  DIO,  400  cps 


Figure  Dll,  46.65  cps 
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Figure  D12,  77.70  cps 


Figure  D13,  88.04  cps 
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Figure  D17,  144.85  cps 


Figure  D19,  188.73  cps 


Figure  D20,  200.05  cps 


Lgure  D21 ,  214.32  cps 


Lgure  D22,  224.78  cps 


Figure  D23,  252.46  cps 


Figure  D24,  259.82  cps 


Lgure  D26,  309.24  cps 


Figure  D28,  330.26  cps 


Lgure  D29,  345.66  cps 


Figure  D30,  372.63  cps 


Figure  D31,  388.01  cps 
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Figure  D32 ,  406.51  cps 
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MODELS  USED  FOR  STEADY 


AERODYNAMIC  ANALYSIS 


NASTRAN  STEADY 
AERODYNAMIC  MODEL 
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APPENDIX  F 

EXPERIMENTAL  AND  NUMERICALLY 
PREDICTED  STEADY  CHORDWISE 
PRESSURE  DISTRIBUTIONS 
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Figures  F1-F6  compare  experimental  and  numerically  predicted 
pressure  distributions  on  the  stab  upper  surface. 


□  -  Experimental 
O-  USSAERO 

F1-F3  a  =  2.3°  M  =  0.19 

F4-F6  a  =  4.8°  M  =  0.19 

Figures  F7-F12  compare  numerically  predicted  pressure 
differences. 

□  -  USSAERO 
O-  NASTRAN 

F7-F9  a  =  2.3°  M  =  0.19 

F10-F12  a  =  4.8°  M  =  0.19 
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This  thesis  concerns  the  development  of  a  finite  element 
model  of  the  T-38  horizontal  stabilator  Cor  use  on  NASTRAN. 

The  model  is  to  be  used  to  analyse  degradations  in  flutter 
speed  due  to  repair. 

Static  analysis  has  shown  the  model  to  be  lacking  in 
torsional,  stiffness.  The  probable  cause  being  the  inability 
of  NASTRAN  plate  bending  elements  to  model  torsion  cells. 
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An  increase  of  elastic  and  shear  moduli  of  plate  bending 
elements  in  the  model  by  30  percent  produced  more  accurate 
results  but  additional  investigation  is  necessary. 

Modal  analysis  has  pointed  to  a  modeling  error  in  the 
root,  trailing  edge  area.  The  affect  has  caused  an  additional 
node  to  appear  on  the  trailing  edge  for  modes  above  100  cps  in 
a  free-free  condition.  Investigation  of  the  steady  aerodynamic 
pressure  distribution  over  the  stabilator  shows  good  correlation 
with  experimental  resuLts. 

A  flutter  analysis  procedure  was  established  and  the 
affects  of  the  errors  found  in  the  structural  model  were 
investigated.  With  no  corrections  made  to  the  model,  a  flutter 
speed  equivalent  to  that  predicted  using  strip  theory  was 
achieved  for  the  sea  level  condition. 
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